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We present a novel method to perform quantum state tomography for many-particle systems which are partic- 
ularly suitable for estimating states in lattice systems such as of ultra-cold atoms in optical lattices. We show that 
the need for measuring a tomographically complete set of observables can be overcome by letting the state evolve 
under some suitably chosen random circuits followed by the measurement of a single observable. We generalize 
known results about the approximation of unitary 2-designs, i.e., certain classes of random unitary matrices, 
by random quantum circuits and connect our findings to the theory of quantum compressed sensing. We show 
that for ultra-cold atoms in optical lattices established experimental techniques like optical super-lattices, laser 
speckles, and time-of-flight measurements are sufficient to perform fully certified, assumption-free tomography. 
This is possible without the need of addressing single sites in any step of the procedure. Combining our ap- 
proach with tensor network methods - in particular the theory of matrix-product states - we identify situations 
where the effort of reconstruction is even constant in the number of lattice sites, allowing in principle to perform 
tomography on large-scale systems readily available in present experiments. 



I. INTRODUCTION 

Quantum state tomography is - for obvious reasons - a pro- 
cedure of great importance in a large number of experiments 
involving quantum systems: It amounts to reconstructing an 
unknown quantum state entirely based on experimental data. 
In many situations one indeed aims at identifying what state 
has actually been prepared in an experiment. This seems par- 
ticularly important in the context of quantum information sci- 
ence, where quantum state and process tomography is now 
routinely applied to small, precisely controlled quantum sys- 
tems (Bill. Yet, needless to say, in a number of other contexts 
the reliable reconstruction of quantum states is an important 
aim as well. 

For finite-dimensional quantum systems, conventional 
quantum state tomography can be performed by choosing 
a suitable basis of S((D^), i.e., the operators on the d- 
dimensional Hilbert space of the system in question. Then, the 
expectation values of these d'^ observables are being measured 
to some required accuracy, from which one can reconstruct the 
unknown density matrix p. The same approach, however, is 
doomed to failure when applied to quantum many-body sys- 
tems: If one has a many -body system at hand with k lattice 
sites of local dimension di, the number of necessary differ- 
ent measurement settings is given by m = df^, i.e., it scales 
exponentially with the size, rendering the treatment even of 
reasonably large systems impossible. Techniques of quantum 
compressed sensing (IHEl allow to significantly reduce the re- 
quired number of measurement settings, if the state is of rank 
r, to m = &{rdlog^ d) (where 6 denotes asymptotic equal- 
ity). If r <C (i, this is an impressive reduction, and gives rise to 
feasible quantum state tomography in medium-sized quantum 
systems, but this number is still exponential in the number of 
sites. Such a scaling cannot be overcome without further re- 
striction of the class of possible states, simply because even a 
pure state needs of the order of d parameters to be described. 
However, if the state is not only pure but also described by a 
generic matrix product state (MPS), the necessary number of 




FIG. 1 . Realization of a random quantum circuit by an optical su- 
perlattice. The lattice with the larger period is switched between the 
two depicted situations, and the lattice depth is changed locally by 
Ai which is uncorrected between the lattice sites. 



measurements only scales linearly with the system size and is 
even independent of it for the important special case of trans- 
lationally invariant MPS. What is more, in several instances 
the classical procedure to reconstruct the MPS matrices from 
the measurement data is efficient 1 9 1 . 

This small number of parameters ought to make tomogra- 
phy an easier task, but in many practical settings involving 
quantum many-body systems, a serious challenge arises: In 
most interesting systems it is very difficult if not impossi- 
ble to directly measure a full operator basis. Instead, merely 
measurements of some preferred observables might be readily 
available. 

In the present work, we propose a solution to this problem. 
We do so by combining the action of a suitable random quan- 
tum circuit with a measurement of a very small number of 
different observables (even a single one can be enough). Such 
random circuits are just becoming a tool of great theoretical 
importance in several subfields of quantum information the- 
ory CMl. 
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Here, we show that they also offer significant technological 
advantages, in that they allow for the natural implementation 
of complete quantum state tomography in systems of ultra- 
cold atoms using only the techniques of super-lattices and 
laser speckles. Both techniques have already been experimen- 
tally proven to be feasible (TSlIISllISl. Combining these new 
insights with the above mentioned methods of MPS tomogra- 
phy brings, for the first time, full tomography of many-particle 
systems close to experimental reach. Note that once the state 
is estimated, one can from this knowledge of course also com- 
pute properties such as entanglement entropies of systems of 
cold atoms 1 14] - but also all other properties that are a func- 
tion of the quantum state. 

Needless to say, it still amounts to a very difficult pre- 
scription. But while conventional measurements in ultra- 
cold bosonic systems amount to measuring certain correlation 
functions or estimates of the temperature in thermometry, say, 
the path described here eventually allows for the full recon- 
struction of an unknown state of a quantum many-body sys- 
tem - a very promising perspective. 

The remainder of this article is organized as follows: First, 
we introduce in Section [Il| tomographically complete sets of 
observables and a generalization of operator bases, called tight 
frames, and discuss their realization by means of random uni- 



tary matrices. In a next step, we show in Section III how ef- 
ficient compressed sensing is possible with families of uni- 
tary matrices which form approximate unitary 2-designs 1 10 - 

fn 



fT7HT9ll . In Section |IVj a way to efficiently realize such 
unitary 2-designs with the help of random quantum circuits is 
presented, before we discuss the application of this approach 
to ultra-cold atoms in optical lattices in Section [V| Finally, 
Section|Vl] shows how, under certain assumptions on the state, 
tomography is possible with a number of measurements which 
is linear or constant in the number of lattice sites before a con- 



clusion is given in Section VII 



II. 



TOMOGRAPHY BY MEANS OF UNITARY 
EVOLUTION 



Before we turn to the question of how tomography is pos- 
sible, we provide general conditions for sets of observables to 
be suitable for reconstructing quantum states. We discuss why 
tomographic completeness as such is not sufficient and intro- 
duce generalizations of operator bases which allow for robust 
tomography. 



operator norm respectively. Furthermore, we use the Hilbert- 
Schmidt scalar product defined as (A, B) = Tt{A^B) and the 
projection on the subspace spanned by some w e S which is 
defined as 



(1) 



Noting that a measurement of the expectation value of w cor- 
responds to determining {w^ p), we define the total sampling 
operator 

Wd = d^ I dp{w)V^ e B{B{n)), (2) 

where /i is a probability measure on S. 

A finite set of observables is said to be tomographically 
complete if any two different quantum states have distinct ex- 
pectation values for some observable: This implies that one 
can theoretically reconstruct the state if one knows all expec- 
tation values. If Wd has full rank, i.e. rank(>V(i) = (i^, the 
state p can be obtained from Wd{p) by matrix inversion if 
issues of statistical errors and numerical imprecisions are ne- 
glected. For general probability measures we make use of the 
subsequent definition. 

Definition 1 (Tomographic completeness). The measure p is 
said to be tomographically complete if Wd, cis defined in Eq. 
(|2]), is full rank. 

In practice, not every tomographically complete measure 
on the observables is necessarily useful: The observables, 
viewed as vectors in , should not be too unevenly dis- 
tributed over the sphere. That is to say, if the ratio between 
the largest and smallest eigenvalues of Wd is large, small er- 
rors in the expectation values can lead to large errors in the 
reconstructed state. The ideal situation is that of a tight frame, 
also known as spherical 1 -design: 

Definition 2 (Tight frame). A probability measure p on the set 
of2-norm normalized Hermitian matrices S is called a tight 
frame ifWd = 1 ^ith Wd given by Eq. (|2]). 

Examples for tight frames include any operator basis and 
the rotationally invariant measure on the 2-norm sphere S . 
When observables are taken from tight frames, the reconstruc- 
tion problem is well conditioned and small errors in the ex- 
pectation values only lead to a small error in the reconstructed 
state [|8il20i. 



A. Tomographically complete sets of observables and tight 
frames 

We consider a quantum system consisting of k subsystems, 
called (lattice) sites which all have local dimension di . Let p G 
5(H) be a quantum state on the Hilbert space H = (C^ with 
dimension d = df and S the set of corresponding normalized 
observables, i.e., Hermitian matrices w G B{l-L) with ||i^||2 = 
1. In this whole work, we denote by || • \\p the Schatten p- 
norm, where p = 2 and p = oo are the Frobenius norm and 



B. Evolution of observables 

A tomographically complete set of observables must con- 
tain at least d'^ observables, which might be difficult to mea- 
sure directly. We introduce a way to do tomography by per- 
forming suitably random unitaries, followed by the estimation 
of the expectation value for a single observable. A different 
but related approach has been employed in Ref. |21 1 to per- 
form high-fidelity quantum state reconstruction in situations 
where the knowledge about the state is not tomographically 
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complete. In the present work, the time evolution is a tool to 
obtain knowledge about the quantum state of the system leav- 
ing the question of determining Hamiltonians aside 1221 . 

Switching to the Heisenberg picture, the outlined procedure 
amounts to measuring the time-evolved observables. Simple 
as this idea is, it allows the economical reconstruction of un- 
known quantum states, as it turns out: 

Definition 3 (Induced observables). To a measure ft on the 
special unitary group SU {d) and an observable wq we asso- 
ciate the following induced measure on S: 



1 



1 



(3) 



where f : U{d) ^ S is defined by f{U) = WwoU, and 
Sx denotes the Dirac in x. Notice that if jl is a probability 
measure, then so is ji. 

An important insight is provided by the following observa- 
tion. 

Tiieorem 1 (Tight frame induced by random unitary). Let w 

be a traceless, normalized observable and /iH be the Haar 
measure on SU{d). The measure induced on S is a tight 
frame. 

Proof: To show this observation, we use ^ to calculate the 
sampling operator ^ for a state p with Tr(p) = 1: 

m^p) = {d^ - 1)/ dMU){U^woU,p)U^woU^^. 

(4) 

As w is diagonalized by a unitary matrix, we can assume it to 
be diagonal and obtain 



X Wk,kU)m,mPl,n 



(5) 



where U ^ pu indicates that U is distributed according to 
the Haar measure. The occurring expectation values can be 
obtained from Ref. |23|. We first consider the off-diagonal 
elements, i.e., the ones with i ^ j. Here, the expectation 
value vanishes unless / = i and n = j, in which case we get 

^^J—^{6k,Un-l) + {l-5k,„,){-l)). (6) 

We now turn to the diagonal elements and note that if i = j, 
one needs / = n to get a non- vanishing expectation value. We 
consider two cases separately: If m = /c we get 

^t/~MH [Uk,iUk,iUk,iUk,i] = 

^^J—^{5,^a{n-l) + {l-5i,i){n-l)), (7) 



while iorm^k we obtain 

Ef/^^H \Uk,lUrn,iUk,lUrn,i\ = 

iki{n-l) + {l-5i,i)n). (8) 

Inserting now ([6]), Jt]), and ^ into ^ and using ^ • pi^i = 1, 
^. Wi^i 0, and f^. |i^i,ir = 1, we get Wd,u{p) P which 
concludes the proof. ■ 



C. Tight frames under physical restrictions 

In many situations of interest, the quantum state is not com- 
pletely arbitrary but satisfies some additional conditions. In 
experiments with ultra-cold atoms, for example, the total par- 
ticle number is conserved and super- selection rules forbid su- 
perpositions of states belonging to different eigenvalues of the 
corresponding operator. Assume that the quantum state p acts 
only in some subspace and denote the corresponding projec- 
tion super-operator as Vn- An important example is given by 
eigenspaces of an operator TV describing a conserved quan- 
tity, i.e. one with commuting with p. In this case of a quantum 
state confined to a subspace. Definition [2] can be relaxed to 



N 



(9) 



where dN denotes the dimension of the subspace correspond- 
ing to the eigenvalue N of N. If one aims at realizing such 
a restricted tight frame by the means of Theorem [T] one can 
replace the group SU{d) hy 

SUj^id) = {[/ G SU{d) : [U,N] = o} . (10) 

As this group is compact, there exists a unique Haar probabil- 
ity measure on it. We adapt Definition [3] to this situation. 

Definition 4 (Induced observable on subspace). Let N e 

B(l-L) and Sn its d^ dimensional eigenspace corresponding 
to the eigenvalue N. To a measure p on the group SUj^{d) 
and an observable wq £ Sn with ||i^o||2 = 1, we associate 
the following measure on Sn' 



1 



1 



1 



^1 



(11) 



where f{U) = U^w^U. 

The matrices p, wq, and all V G SUj^{d) are block di- 
agonal with respect to the eigenbasis of TV. We consider the 
block corresponding to the eigenvalue N. Invariance of the 
Haar measure on SUj^ implies that if V ~ Mh,^^/^^^), then 
Vn{V) ~ P}i,su{dN)- Thus, one can apply the proof of The- 
orem [T] to this block and obtain Eq. 



III. UNITARY t-DESIGNS 

The effort to implement random unitaries drawn from the 
Haar measure scales exponentially in the number of lattice 
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sites k, making an implementation both theoretically ineffi- 
cient and pratically unfeasible. However, this problem can be 
circumvented by replacing the Haar measure by a unitary 2- 
design which is much easier to sample from as we will see 
later on. Unitary t-designs behave like the Haar measure in 
specific situations |[l2|[T8||T9l. The definition most suited to 
our problem is the following: Let uhc 3. probability measure 
on SU{d). We define two channels on B{n^^): 



(12) 



and Qt^nn^ where /in is the Haar measure. We say that u 
is a unitary t-design if Q^, = Qn. We say that v is an s- 
approximate t-design if 



where || • | 
defined as 



WGt^u - Gt^uW < e (13) 
denotes the superoperator 2 2-norm which is 



\0\\ = sup 

^,11^112 = 



\\o{x)h 



(14) 



This superoperator norm is equal to the Schatten oo-norm 
when the channel is seen as a mere linear operator acting on 
the real vector space of Hermitian matrices. In the remain- 
der we only consider the case t = 2 and drop the index t for 
simplicity. 



A. Tight frames from unitary 2-designs 

For our purpose, i.e., replacing the Haar measure random 
unitary in Theorem [T] we need an approximate unitary 2- 
design. This is the case because when Haar measure induces 
a tight frame, c.f. Eq. both U and W appear twice. 

Theorem 2 (Tight frames induced by unitary 2-designs). Let 

Wo be a traceless, 2-norm normalized observable and v be an 
e-approximate unitary 2-design. Then the sampling operator 
Wd,u corresponding to the induced measure fulfills 



\\Wd,u-n<Vd{d^ -l)e. 



(15) 



(16) 



Proof: We have 

= sup ||Wrf,.(X)-W^,H(X)||2. 

^,11^112 = 1 

We note that 



mA^) ={d^ - l)Et/^. [{U^woU,X)U^woU] + -TrX 




TrX, 



(17) 



where Tri denotes the partial trace with respect to the first of 
the two subsystems of equal dimension. Note that we have ex- 
tended the definition of yVd,u to operators with non-unit trace. 
This relation yields ( p3) ) after inserting it into Eq. (16) and 
applying ([13]). ■ 
The same argument holds also for restricted tight frames as 
defined in Let be a distribution on SUjy{d) such that if 
U ^ V, then Vn(U) is drawn from an e-approximate 2-design 
on SU{dN)- This implies 

Wm^^uVN -Vn\\< Vd^{d% - l)e, (18) 

as follows from the application of the above proof to the block 
in U corresponding to the eigenvalue N of N. 



B. Compressed sensing 

The technique of compressed sensing allows to reduce the 
number of measurements which are needed to reconstruct a 
quantum state from Q{d'^) to 9((ipolylog((i)) if the rank of 
the state does not increase with d. To perform this method, one 
has to choose m = 6((ipolylog((i)) observables . . . , Wm 
randomly from the tight frame according to the correspond- 
ing probability measure and determine their expectation value 
{wi^ p) by measurement. Then, one can efficiently solve the 
optimization problem 

min ||cr||i s.t. Vz = 1, . . . ,m : {wi.cf) = {wi, p). (19) 

(7 

The theory has been developed for observables forming oper- 
ator bases in Refs. (SlIU and extended to tight frames in Ref. 
|8|. There, it was also shown that the tight frame condition 
may be violated and compressed sensing is still possible if 



111 < 



1 



8v^' 



(20) 



where r is the rank of the state we want to reconstruct. Not 
all tight frames are equally suited for compressed sensing as 
can be seen with a simple example: Let . . . , k;^2 be the 
elements of an orthonomal operator bases of S(C^) where wi 
is a rank-one projector and p = wi. In this case, one has to 
measure of the order of d'^ observables before one "hits" wi 
and gets any information on the system. In Refs. |[5l|6l[8l, it 
has been shown that this problem cannot occur if all observ- 
ables fulfill the so-called "Fourier type incoherence condition" 
which reads 



P. 



I^ll 



A 



d 



0, 



(21) 



where A must fulfill A = O (poly log (d)). Note that statements 
like ( [2T] ) make only sense when considering families of tight 
frames with growing dimension d. As we are mainly inter- 
ested in the asymptotic efficiency of our scheme, we restrict 
ourselves to the scaling behavior and omit explicit pref actors. 
We now give a condition under which Eq. (21 ) is fulfilled in 
the situation of interest. 
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Theorem 3 (Compressed sensing with induced observables). 

Let Wo be a traceless, normalized observable fulfilling 
ll'^oll^ ^ ^/d with A = 0(polylog((i)), and let v be a 
l/{SVrd{d'^ — 1)) -approximate 2-design. The induced tight 
frame fulfills ( [27] ), which implies that it allows for compressed 
sensing. 



Proof: Since 



(22) 



and ||1/V^||^ = l/d, condition (21) is fulfilled and from 
( p3] ) it follows that ( [20| ) is satisfied which proves that com- 
pressed sensing is possible. ■ 
Theorem [3] holds especially in the important case of a ob- 
servable which acts non-trivially only on a few number of lat- 
tice sites because here it is of the form 



1 ik — m 

(4-^)1/2 



(23) 



for some normalized, traceless observable v and some small 
constant m. 



IV. APPROXIMATION OF UNITARY 2-DESIGNS BY 
RANDOM QUANTUM CIRCUITS 

In the previous section, we have shown that unitary 2- 
designs can be used to realize tight frames. We now show how 
they can be approximated by parallel random circuits and gen- 
eralize the results of Refs. ifTOlfTTlfTTl to show the following: 

Theorem 4 (Random circuits). Assume k to be even. Con- 
sider a parallel random circuit where in each step either 

^1,2 ^ %,4 (8) ... Uk-l,k or 1/2,3 ^ ^4,5 ... (8) Uk-2,k-l 

is performed with probability 1/2 where Ui^i-^i acts on the 
neighboring sites i and i -\- 1. If the nearest-neighbor uni- 
taries are drawn, in each step independently, from a probabil- 
ity measure V2 which is universal, as defined below, there ex- 
ists a constant C (depending on the local dimension) such that 
the random circuit forms an e-approximate unitary 2-design 
after n = C \og{l/e)k log k steps. 

A finite set of nearest-neighbor unitary quantum gates is 
called universal if they generate a dense subgroup of SU {df). 
For an arbitrary probability measure on SU{df), the notion of 
universality can be generalized, according to Ref . fTl\ : 

Definition 5 (Universality). We say that ja is universal if for 
any open ball S there exists I > such that S has a nonzero 
weight for the l-fold convolution product of fi. 

Proof: The proof of Theorem |4] is an extension of that of 
similar results in Refs. ifTOlfTTl . Readers mostly interested in 
the application to optical lattice systems can safely skip to the 
next section. 

In Ref. fTPl, it is shown that parallel random circuits with 
periodic boundary conditions form £-approximate 2-designs 
after n = C\og{l/e)k steps if the unitaries are drawn from 
the Haar measure on SU{df). 



We now proceed in three steps: First, we show that the 
nearest-neighbor unitaries can be drawn from an approximate 
2-design on SU (df) instead from the Haar measure. We de- 
note the measure corresponding to a single step of the random 
circuit by and define the linear operator G^^^ by 



This operator can be decomposed as G^,, 
with 



(24) 



(Me + Mo)/2 



Me = Pl,2 (8) P3,4 (8) . . . (8) P/c-l,/c, 
Mo = P/c,l P2,3 8) . . . (8) Pfe-2,fe-l 



where 



P 



We have to bound 



\\GZk - GhIIoo = >^2[Gy^Y 



(25) 
(26) 



(27) 



(28) 



where n is the depth of the circuit and A2 denotes the second 
largest eigenvalue. In Ref. 1 1 1 1, it is shown that if the nearest- 
neighbor unitaries are drawn from the Haar measure, there 
exists a constant A > such that the corresponding operator 
Gy^ fulfills \2{Gy^) < 1 — A. Using now the fact, that ^2 is 
a (^-approximate 2-design, we get 

A2(aj < A2(aj + fe||G,,-GHj|oo < 1 + M-A. (29) 

For the right-hand side to be smaller than one, which is nec- 
essary and sufficient for an exponentially fast convergence, 
one has to choose 5 = 0{l/k). To realize this local approx- 
imate 2-design, we use a result from Ref. [|17il which states 
that one needs to draw only s = 0(log(l/^)) gates from an 
arbitrary universal gate set to achieve this. We pick 5 to be a 
power of 2, which is surely always possible. Thus, we have a 
circuit with depth 



= 0{ns) = 0{\og{l/e)k\ogk) 



(30) 



where the random choice between Mg and Mo is not made in 
every step but in blocks of s steps which corresponds to the 
operator (M| + M^)/2 while s steps of the actual quantum 
circuit performed are described by ((Mg + Mo)/ 2)^. They 
both have the same fixed point. What is more. Theorem [5] 
holds true, see below. This means, in particular, that the con- 



vergence of the actual circuit cannot be slower and ( 30 ) holds. 

The last thing needed to obtain Theorem [4] is to switch to 
open boundary conditions, i.e., remove the first tensor factor 
in Eq. ([26]). As this does not affect the fixed point and the op- 
erator norm difference is on the order of S, only the pref actor 
is changed slightly. We note that the prefactor depends on the 
actual choice of 1^2 • 

The conditions for Theoremfsjto apply are fulfilled if 1 /e = 
0((i^/^). Using this in Eq. ([s^we get 



riT = 0{k^ log k). 



(31) 
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This means compressed sensing is possible with a single trace- 
less observable and a parallel random quantum circuit with a 



depth given by Eq. ( 31 ). 

We did not explicitly discuss the case of restricted tight 
frames because it can be, again treated by block-decomposing 
all matrices according to the spectral decomposition of the op- 
erator N describing the symmetry, c.f. Eq. |iq| ). Thus, with 
a parallel random quantum circuit as in Theorem]?] with uni- 
taries which are universal for SUj^{df) with a depth as in Eq. 
( [3T| ), one can perform compressed sensing for states within 
some eigenspace of TV. ■ 

Theorem 5 (Mixing properties of circuits). Let s be a power 
of 2, then 



Proof: We now prove the validity of Eq. ( |32| ) for s being 
a power of 2, i.e. 5 = 2^ . Let A and B be two Hermitian 
matrices. From < {A — B)'^ it follows that 



A^BY ^ A^^B^ 



2 J - 2 
and hence also, from the monotonicity of eigenvalues 1281 . 



(33) 




(34) 



We also need the fact that, for all pairs of positive Hermitian 
matrices A and B, 



\2{B) > X2{A) \2{B') > X2{A'), 



(35) 



which can be easily seen by diagonalizing A and B. Note that 
Mg and are positive operators. For any natural n, ^2 one 
finds that 



2*2' 



<A2 




Mr + Mr ' 



2*2- 



V. OPTICAL LATTICE SYSTEMS 

We now introduce a novel method for quantum state tomog- 
raphy for systems of ultra-cold atoms in optical lattices that 
does not require any local addressing of single sites. In fact, 
each of the steps necessary has already been demonstrated 
experimentally. The idea is to make use of appropriate ran- 
domness and natural time evolution, suitable exploiting opti- 
cal superlattices, such that time-of-flight images give rise to 
complete tomographic information - quite an interesting and 
promising perspective. Using these experimental tools, the 
above mathematical methods become applicable. It should be 
clear, as pointed out before, that such a prescription also al- 
lows for detecting entanglement entropies 1 14, 29] in systems 
of cold atoms. 

Ultra-cold atoms in optical lattices form some of the clean- 
est quantum many -particle systems available for experiments 
and allow for the realization of well-known effects from con- 
densed matter. For example, both the bosonic super-fluid 
Mott-insulator transition and the Mott state of fermions where 
observed by changing the intensity of the laser forming the 
lattice |24ll25|. Such systems also have the potential of func- 
tioning as quantum simulators which means, they allow to 
simulate systems from other branches of physics like the no- 
toriously difficult quantum chromodynamics (QCD) 1.26 J . 



A. Time-of-flight measurements 

Even though measurements with spatial resolution have 
been demonstrated in recent experiments |16|, the standard 
technique is still provided by time-of-flight absorption imag- 
ing. Here, the lattice and the confining trap are instanta- 
neously switched off and, after some time during which the 
atoms expand approximately without interaction, an absorp- 
tion image is taken [24, 30 1. As the distance that the atoms 
fly during the expansion is proportional to the initial momen- 
tum, this procedure amounts to a measurement of the den- 
sity in momentum space. We restrict ourselves to the bosonic 
case, while noting that fermions can be treated in a completely 
analogous way, and expand the field operators of the one- 
dimensional bosonic field as 



2n+i 



9*2"!" 



(38) 



(36) 



where we have employed ( 35 ) and ( 34 ) for A = ^ and 
B = Mf^ . By a repeated application of (36(, starting from 



ii = 1 and i2 = j, until ii = j and i2 = j, we finally get 



(37) 



where W'^^^ is the Wannier function of the j-th band, is the 
position, and bi^-' the corresponding annihilation operator at 
site s. If the lattice is sufficiently deep, all bands but the lowest 
one can be neglected, and we drop the upper index in ( [38] ). 
The momentum- space distribution, which is measured in the 
time-of-flight experiment, is given by n{p) = \W{p)\'^S{p) 
where W is the Fourier transform of the Wannier function and 
the quasi-momentum distribution is given by 



which is the statement to be shown. 



s{p) = E 



Jp(Xs-Xi 



(39) 
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As we do not assume the state to be translationally invariant, 
Eq. ( [39| ) cannot be inverted to get the two-point correlation 
functions in real space but we can get for integer / 



1 

2^ 



(40) 



where we have set the lattice spacing to one. If the atoms are 
bosons, the local Hilbert space is infinite-dimensional. How- 
ever, as the interaction between the atoms must always be re- 
pulsive to ensure stability of the quantum gas, one can neglect 
state where more than a given cut-off number of atoms are 
present on a single lattice site. This allows us to work with 
a finite-dimensional Hilbert space. We note that this does not 
even needs to be an approximation as one can set the maximal 
number of bosons per site Ns to their total number N. How- 
ever, in any practical setting, one would use Ns <C N and 
still get a very good approximation. For an arbitrary i, we set 



(X 



(41) 



3 = 1 



which is traceless. Due to the sum, which stems from the 
absence of translational invariance, Eq. ( [4T] ) is not exactly of 
the form given by ( [23] ) but a sum of few, i.e. logarithmically 
many in the Hilbert space dimension d, terms of this form. 
This implies 



ik'^iioo<c 



k 

El 



< Ck 



(42) 



where C, C > are constants. Because k = 9(log(i), we 
can employ Theorem [T] to show that a measurement of the 
momentum space distribution, together with an approximate 
2-design, allows for efficient compressed sensing. 

Although already a single choice of i yields an approxi- 
mate tight frame, we can use the data corresponding to all 
i = 1, . . . , /c, as they are measured anyway, to reduce the nec- 
essary number of experiments. 



B. Realization of the random circuit 

We now discuss how a probability measure on the nearest- 
neighbor unitaries that is universal can be obtained. To be as 
specific and simple as possible, we use the single-band Bose 
Hubbard model with Hamiltonian 131 1 



H 



^ U 
E Ji(^l^i+i+^l+i^i) + ^n,(n,-l) + A,n,, (43) 



where hi = blbi and Ji . . . , J/^; Ai, . . . , A^; Ui,...,Uk e 
R. To realize the parallel random quantum circuit, we make 
use of the techniques of super-lattices |[T5]|T6|[27l and speckle 
patterns 1 13 1. As the total number of atoms is conserved and 
super-selection rules forbid the superposition of states with 



different particle numbers, we restrict ourselves to particle- 
number conserving operations. 

By using an additional lattice for which the lattice constant 
is twice as large, one can change the height of the wells be- 
tween alternating pairs of site. This mainly affects the hopping 
constants Ji and to much less extent the interaction parameter 
Ui, an effect which we neglect. By choosing the depth of the 
super-lattice large enough, we get Ji = for the non-coupled 
pairs and Ji = J for the coupled ones. Such a double- well 
structure has been used in Ref. 1 16] to probe correlation func- 
tions. A speckle pattern is created by illuminating an uneven 
surface with a laser and can be modelled by a spatially fluctu- 
ating A^. The situation is sketched in Figure[T] For reasons of 
simplicity, we only consider the regime for which the strength 
of the speckle potential is not correlated over the lattice sites, 
i.e., all A^ are independently distributed. Thus, we have a 
random circuit as in Theorem |4] and we only need to show 
that the corresponding gates generate a dense set in SUN{df). 
The local gates are 

[/,,,(Ai,A2,t) =exp{-it{-J{b\b2^blbi) 



+ ^(ni(ni - 1) + n2(n2 - 1)) 



• Aini + A2n2)), 



(44) 



where t > is the time after which the super-lattice is 
switched and a new realization of the laser speckle is cre- 
ated. We assume Ai 2 to be Gaussian distributed, and we 
have neglected global phases. Now one can adopt an argument 
used in Ref. 1 32 1 for showing that the Gaussian operations to- 
gether with a single non-Gaussian one allow for continuous- 
variable quantum computation. For universality to hold, one 
has to generate all operations where the corresponding Hamil- 
tonian is a polynomial in the creation and annihilation opera- 
tors where every monomial must contain an equal number of 
creation and annihilation operators, i.e., must be balanced, to 



ensure particle-number conservation. This is true as Eq. (44) 
contains all quadratic terms and a single quartic one. Since 
one can generate the entire algebra generated from the origi- 
nal set of Hamiltonian by commutation, one can approximate 
an arbitrary unitary 1321 . Thus, by varying Ai^2, we can ap- 
proximate any gate to arbitrary accuracy which implies, by 
continuity of Eq. (44) universality as in Definition [5] By ap- 
propriately choosing the distribution from which t is chosen, 
the set can be made closed under Hermitian conjugation. Now, 
we can apply Theorem [4] which shows together with Theorem 
[T] that one can perform efficient compressed sensing by us- 
ing a optical super-lattices, laser speckles, and time-of-flight 
imaging. Again, no single site addressing is necessary, and 
still complete tomographic information is obtained. 



VI. MORE EFFICIENT TOMOGRAPHY SCHEME FOR 
MATRIX PRODUCT STATES AND OPERATORS 

Even though compressed sensing notably reduces the num- 
ber of necessary measurements, it still scales exponentially 
with the number of lattice sites. Without any further assump- 
tion on the state, this cannot be overcome. However, when 
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the state is described by a generic matrix-product state (MPS) 
with a fixed bond dimension, tomography is possible with the 
number of measurements growing, in general, almost linearly 
with the system size and being constant for translationally in- 
variant MPS. Ground states of gapped Hamiltonians which 
are a sum of terms acting only on a constant number of lattice 
sites are, generically, of this type |[33ll . 



A. Reconstructing reduced density matrices 

The exponential reduction of the necessary number of mea- 
surements if the state is a MPS is due to the fact that such 
states are completely determined by their reduced density ma- 
trices on all blocks of / consecutive sites where / only depends 
on the bond dimension |i9i|33j|. In Ref. |9|, an efficient algo- 
rithm is given for finding the MPS matrices from these re- 
duced density matrices. Note that this procedure only works 
if an upper bound to the bond-dimension, or, equivalently, to 
the locality size of the Hamiltonian is known. This is obvious 
as there could always be long range correlations which do not 
affect the / site reduced density matrices. 

Let dB = d\he the dimension of the subsystem under con- 
sideration and define the operator Tr^ : S(C'^) ^ S(C'^^) 
be the operator acting as Tr^ (p) = T^^r^P and Rq denote all 
lattice sites but . . . , g + / — 1. As we only want to perform 
tomography on the lattice sites + / — 1, we can trace 

out the remainder of the system. In this case, the tight-frame 
condition of Definition |2] becomes 

TR^oWd,=TR^. (45) 
If additional constraints apply, the condition reads 

Tr^ o Wd^^ o Vn = Tr^ o Vn (46) 

where dsN is the dimension of the matrix block correspond- 
ing to the eigenvalue of TV, restricted to the subsystem of / 
lattice sites. We concentrate on the former case as the second 
follows in an analogous way and show that a random circuit 
with a depth which does not depend on the system size can 
realize such a reduced tight frame: Assume wq to be a sum of 
traceless observables, each acting non-trivially only on some 
block of / lattice sites. For compressed sensing to be possi- 
ble, it is sufficient for v to induce, for every an approximate 
reduced tight frame with 



l)e 



(47) 



where e must be chosen such that y^~dB{(fg — l)s < l/(8\/r, 



c.f. Eq. ( (20j i. This is the case if it is an approximate reduced 
unitary 2-design with 



TTTTTT 



FIG. 2. Top: Parallel random circuit acting on the entire system. 
Bottom: Random circuit used in the proof of Theorem [6] which acts 
only on the inner sites. The above circuit randomizes the states of 
the inner sites not less than the above one. 



To see that this is true, we calculate 



sup 

^,11^112 = 



||Tr^,(W,,,.(X) - W,,,h(X))||2, (49) 



which yields the desired result after inserting (17 ) and apply 



ing(|48). 

To show that obtaining such a tight frame is efficiently pos- 
sible, we adapt Theorem |4j 

Theorem 6 (Reduced tight frames by reduced unitary 
2-designs). Let the parallel random circuit be as in Theo- 
rem^ There exists some constant C such that it forms an 
e-approximate reduced l-site 2-design as defined in ( [?7| ) after 
n = Clog{l/e)l log / steps. 

Proof: We show the fast convergence of our random cir- 
cuit to a reduced unitary 2-design by comparing it with an- 
other random circuit which is easier to deal with, see Fig. |2] 



In analogy to Eq. (|24|) we denote by Olf the channel corre 



sponding to an application of the parallel quantum circuit to a 
block consisting of / lattice sites starting from q. Theorem |4] 
implies 



sup WTTR^iiGify 

X,\\X\\2 = 1 



(50) 



for n = Clog{l/e)llogl where C is a constant. Using that 
^2{Guk) ^ ^2{Gui), one obtains (48), concluding the proof. 



: sup 

^,11^112 = 1 



TTR^{g.-gl^^){x) 



< e. 



(48) 



where denotes the channel which acts as on a block of 
/ sites starting at q and as the identity on the rest of the system. 



B. Complexity of classical post-processing 

Theorem [6] implies that a random quantum circuit of con- 
stant depth is sufficient to perform tomography on a reduced 
density matrix of constant size. The number of operations, i.e. 
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FIG. 3. Influence region for a local observables. Only the action of 
the darkly colored gates influences the measured observables. 



random unitary gates and measurements of expectation val- 
ues of wq, does only scale polynomially with the number of 
lattice sites. Therefore, we regard the quantum part of the 
protocol as efficient. However, this says nothing about the 
amount of post-processing needed because the reconstruction 
of the /-site reduced density matrices requires the knowledge 
of Ttr^ w for all observables w obtained by the realizations of 
the random quantum circuit. If one needed to keep track of the 
evolution of observables on the entire Hilbert space this would 
require an exponential amount of computational resources. To 
see that this is not a problem in the present situation, we use 
the fact that wq is a sum of terms which act non-trivially only 
on blocks of constant size. Thus, the observables induced by 
the constant-depth random circuit can be written as 

w{U) = Y,w,{U) iRjidR^^^ (51) 

i 

where Wi{U) acts non-trivially only on a block of constant 
size L starting with site i and where Ri denote all lattice sites 
but i, . . . , i + L — 1, see FigurejS] From Eq. ( [5T] ), we get 

Trfi,(w;(/7)) = ^ v^Trfi,K(f/)), (52) 

which means that one only needs to deal with observables on k 
Hilbert spaces which all have dimension dj^ not depending on 
k making also the classical part of the protocol efficient. If the 
system is assume to be translationally invariant, all reduced 
density-matrices are equal, reducing the necessary number of 
measurements to a constant. Roughly speaking, the random 
circuit transforms a local observable to a reduced tight frame 
on / lattice sites with some influence on L > / sites and none 
on the rest of the system. 

C. Mixed states 

Even though the method developed in this section is, in its 
present form, limited to pure states, it can be naturally ex- 
tented to mixed states by using recent results on tomography 



for matrix-product operators (MPO) which are a natural gen- 
eralization of MPS [34.. 35 1 . In Ref. [36J, it is shown that large 
classes of MPO states can be efficiently reconstructed from 
the reduced density matrices on a constant number of sites 
not depending on the system size. As we have presented a 
method of recovering these objects, one directly obtains a way 
of performing tomography on states in optical lattices which 
are described be MPO while requiring the same experimental 
techniques. 



VII. CONCLUSION 

In this article, we have presented a new route towards effi- 
cient quantum state tomography for quantum many-body sys- 
tems, specifically for bosons in optical lattices. By using ran- 
dom circuits, which can be implemented by means of super- 
lattices and laser speckles, one can avoid the use of tomo- 
graphically complete local measurements on single sites (in 
fact any local addressing) - still complete tomographic knowl- 
edge can be achieved. These are challenging, needless to say, 
but rely solely on time-of-flight imaging techniques, which 
are nowadays routinely implemented. Without any further as- 
sumptions to the state, the number of necessary measurements 
is optimal up to constants and logarithmic factors in the sys- 
tems dimension. Restricting the set of possible states to ma- 
trix product states, both the number of measurements and the 
depth of the required random quantum circuit does not at all 
depend on the system size. This idea gives rise to the exciting 
perspective of actually reiably measuring out the full quantum 
state of a quantum many-body system in the laboratory. 

There are a number of questions arising from this: For ex- 
ample, it would be very interesting to compare the perfor- 
mance of the present scheme, which is based on random cir- 
cuits, with one where the quantum gates are chosen in an op- 
timal way from some set of feasible operations. While there 
is not much room for improvements concerning the asymp- 
totic behavior, the performance for small systems might differ 
notably. This will be studied extensively by numerical means 
in forthcoming work. We hope that this work stimulates fur- 
ther work - both of theoretical and especially experimental 
kind - in "quantum system identification", in order to inno- 
vate ways of rendering quantum state tomography feasible for 
large quantum systems. 
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